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Abstract 

Radiative capture of the K~ by the deuteron is examined as a reaction for measure- 
ment of the A-neutron scattering lengths. Using the final state interaction, analogous 
to measurements of the radiative capture of pions by deuterium, the scattering lengths 
of both the triplet and singlet A-n interaction can be inferred. The problem of the 
separation of these two fundamental parameters without and with spin information in 
the experiment is addressed. It is shown that a measurement of the deviation of the 
vector deuteron polarization asymmetry from -1 provides a direct determination of the 
difference of the singlet and triplet scattering lengths. 
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1 Introduction 



The interaction of the nucleon with the A hyperon is of great theoretical interest. Since 
the one-pion-exchange interaction is absent (except under broken isospin |l]]), one is 
able to focus on the shorter range contributions to the interaction. The singlet and 
triplet scattering lengths are predicted to be equal by SU(6) symmetry. 

Unfortunately, experiments are difficult to perform for this system because of the 
short lifetime of the A. Nevertheless, some experimental results do exist Q and have 
been analyzed^, Rijken et al. || have recently generated a set of potentials which 
fit these data. The range of scattering lengths obtained from this analysis is rather 
large, which is perhaps a natural consequence considering that the data do not extend 
to low energies. 

Use of the radiative capture of a 7r~ from an s-wave atomic state by deuterium to 
measure the neutron-neutron (nn) scattering length, through the influence of the final 
state interaction, was investigated in an experiment by Phillips and Crowe Q and soon 
after studied theoretically by McV6y||] and Bander Q. The basic technique consists 
of comparing the shape of the of the photon spectrum with calculations. Because of 
the final-state interaction of the two neutrons, a peak results near the high-energy 
end of the spectrum corresponding to a maximum in the s-wave phase shifts at low 
neutron-neutron energy. Since this phase shift rises very quickly to a maximum, as 
a function of energy in the neutron-neutron frame, the peak in the photon spectrum 
is very near the maximum energy. It is also possible to use the neutron spectrum at 
low energies for the studies provided that the angle between the neutron and photon 
is known. Experiments have been performed in both geometries. 

These early exploratory efforts were followed by an experiment detecting the neutrons []|] 
which gave useful results. Following a study of limits due to the theoretical uncertainties ||] , 
experiments were carried out at the Paul Scherrer Institute (PSI) detecting only the 
photon, as well as experiments in which the photon and neutron were detected in 
coincidence |10[|. These experiments resulted in values for the nn scattering length (and 
effective range) which were considered to be the best obtained to date. Recently an- 
other experiment was performed using this reaction at the Clinton P. Anderson Meson 
Physics Facility [|ll|]. Results from this experiment appear to confirm the PSI results. 
The uncertainty in the scattering length is of the order of ±0.5 fm for a value of -18.58 
fm, about 2-3%. 

A determination of the An scattering lengths from the photon spectrum of the 
reaction K~d^nA7, where the capture takes place from an atomic state, was suggested 
by Gibson et al. [|l^] . A study of this process performed by Workman and Fearing fllf| 
concluded that the different hadronic routes considered by Akhiezer et al. [fbtl had a 
negligible impact on a possible measurement (however, see Refs. for comments 
on their representation of the amplitudes). They also concluded that the dominant 
operator for the conversion of a proton to a A through radiative capture of a kaon was 
of the Kroll-Ruderman form 
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= cre (1) 

where e is the polarization of the outgoing photon and <x is the spin operator for the 
proton on which the capture takes place. 

At the same time as this theoretical effort, a feasibility study [ 16 demonstrated that 



the spectrum for this reaction could be separated from the background. Workman and 
Fearing compared their analysis with these data to show that reasonable values of the 
scattering lengths were consistent with the data. 

Two important issues were left unresolved by the work of Workman and Fearing. 
They used only the asymptotic A- neutron wave functions. In the work of Ref. |J 
the lack of knowledge of the short range final-state wave function was found to be the 
dominant theoretical uncertainty in the analysis of the scattering length, of the order 
of ±0.3 fm. If the uncertainty is of the same order in the case of A-neutron scattering 
the consequences are more serious since the scattering length itself is thought to be an 
order of magnitude smaller. 

A second problem, which may be more serious, is that, while in the case of the 
neutron-neutron final state the spin state of the two neutrons is restricted to be singlet 
by the Pauli principle, the An final state consists of a mixture of a singlet and triplet 
states. One expects that the scattering in these two states may be different, and indeed 
the dependence of the scattering length on the spin might provide some very important 
clues to the structure of the interaction. While it has been suspected that polarization 
information may provide a tool to obtain a singlet-triplet separation (see Balewski et 



al. [17] for example), we provide here the formalism and expressions necessary to obtain 
such a separation for this reaction. 

We first give a brief heuristic description of the basic physics underlying the devel- 
opment to follow. Consider a deuteron target prepared in a configuration with only 
magnetic spin projection zero along a z-axis defined by the direction of the photon. 
Assume also that the final state of the A-n system is in a singlet state. Since the K~ 
is assumed to be captured from an atomic s-state, there will be a total spin projection 
of zero in the initial state and zero angular momentum projection in the final state, 
other than that due to the spin of the 7. Since the "transversality" relation, 



ek = 0, (2) 

requires that the spin projection of the photon must lie along its direction of travel (i.e. 
there is no zero spin projection) this assumed transition to the singlet final state must 
have zero amplitude due to the conservation of the z-projection of angular momentum. 
Thus, for example, a spectrum taken under these conditions could be analyzed to 
obtain the triplet scattering parameters alone. 

The previous paragraph gives only an example of how deuteron spin information 
can affect the analysis. A measurement of the photon circular polarization can also 
help to separate the singlet and triplet scattering states. In the subsequent sections 
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we develop the expressions necessary to calculate the photon spectrum and provide 
specific results. 



2 Expressions for the Rate 

For comparison with the measurement we need the theoretical prediction for the shape 
of the spectrum. The absolute rate is not usually measured in this type of experiment. 



2.1 Kinematics and Phase Space 

The magnitude of the momentum, p, of each baryon in this pair's center of mass is 
given by 

pl J, + Ml-Mlf _ Ml (3) 

with s the center of mass energy of the A-neutron pair given by 

s = E%-2kE , (4) 

where Eq is the total energy of the initial system and k is the final photon momentum. 
M\ and M n are the masses of the A and neutron respectively. 

Since we shall assume that the final-state baryons are not observed, we must inte- 
grate over their directions, giving for the differential rate of the observed photon 

dT oc^^/|M| 2 ^ p (5) 



dVt^dk (uj\ + uj n ) 

The total energies, tu n and lv\ are very nearly constant over the energy range of interest. 

It should be noted that a single factor of the photon momentum k appears. It is 
linear only with relativistic phase space; it would appear to the second power in non- 
relativistic phase space. When this reaction was first proposed, the non-relativistic 
version was usually used (in fact Williams |lfi] used this reaction as an example of the 



use of non-relativistic phase space). Of course, either type of phase space can be used 
provided that the proper corresponding operator is used. The difference between a first 
and second power of the photon momentum is relativistic in origin and the effect is 



very small in many cases. In the analysis of Ref. [10] only the very end of the spectrum 
was used and the photon momentum is sufficiently constant over this range that such 
a factor is immaterial. For the analysis of Ref. jnj the larger range of the spectrum 
used leads to a small sensitivity to the power of k used in this expression. For the 
present reaction, since the peak lies farther from the end point of the spectrum than 
for the ir~d — ► nnj reaction, a larger range of photon momenta may be included in 
the analysis and the correct factor is more important. 
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2.2 Matrix Element 



The mathematical development of the expressions for the matrix element used here 



follows a different procedure than that of Refs. [|| and 13]. The previous methods 
have written the matrix element for the transition first as a plane wave transition 
then corrected this expression by adding and subtracting the s-wave contribution. The 
s-wave contribution with a final state interaction is then substituted in the added 
element. For a matrix element with two particles detected in the final state this is the 
appropriate method. However, for the case in which only the photon is detected, there 
is an average over the unobserved baryon momenta which must be made as outlined 
in the previous section. This integral was carried out numerically in the previous 
developments. 

In the method used in the present analysis the partial-wave quantities are calculated 
in terms of the relative A-neutron momentum. This technique requires several partial 
waves to be computed for a handful of magnetic quantum numbers. While it may 
seem a priori less efficient, in the end it has several advantages. The first is that the 
final averaging over the direction of the unobserved momenta can be done analytically. 
Since the contributions from the higher partial waves depend only on the initial state 
wave functions (we assume no interaction in partial waves higher than £ = 0), if one 
wishes to fit the data by varying the scattering length and effective range for the singlet 
and triplet states, these higher partial waves need be calculated only once. 

A second advantage is that the degree of coherence or incoherence of the singlet 
and triplet states is manifest. For example, we will see that if all magnetic quantum 
numbers of the initial deuteron are equally populated and the polarization of the fi- 
nal photon is not observed, the transitions to the singlet and triplet final states are 
incoherent. 

The matrix element for the reaction which proceeds from an initial deuteron with 
spin-projection S z to a final state of the A-neutron system in a total spin state (S' , S' z ) 
is 



jyj-S' ,S' Z ,S Z 



J dr 1( fr 2 < 5 / ^|^ 5 ,(p,r)e- iK -( ri+r2 )|Oe ik - ri |^(r) > (6) 



where ipg, (p, r) is the final state wave function of the neutron- A system and the variable 
S' is or 1 for the singlet or triplet case. K is the center-of-mass momentum of the 
baryon pair. The quantity |^^f(r) > is the deuteron wave function with an initial 
spin projection S z . We obtain the radial deuteron wave function from the solution 
of the coupled Schrodinger equations for a pure one-pion-exchange potential. It has 
been shown that this interaction reproduces all of the low-energy properties of the 
deuteron [19, 20, ^] (likely to be the most important in this calculation). The operator 



= o-i e. 

Transforming to relative and center-of-mass coordinates in the spatial variables we 
have 

M S',S' Z ,S Z = ( 27 r) 3 (5(K - k) J dr < S^|^,(p, r)|0e^ |tf£f (r) > . (7) 
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We write the deuteron wave function as 

|*«T(r) >= S(r)Y °\lS z > +D(r) £ Y 2 m ' (r)\la > (8) 

and the A-n relative motion wave function as 

V>5'( P ,r) = 4nJ2i% m (r)Y e m * (p)*f(r). (9) 

With these expansions we can write (with an arbitrary overall normalization, omit- 
ting the delta function) 



M s '^=4n<S'S' z \0\ [dr £ i^V* (r)Y L M (r)Yr(p)Y L M *(k)4>f (r)j L (% 

•> o ™ r »,f V Z 



S(r)Y °|lS 2 > +J D(r)^y 2 m '(r)|la > 



With the photon direction, k, along the z axis, one has 



y L M (k) = 5 mm 



'2L + 1 



47T 



(10) 



(11) 



Performing the integral over the angles of r, we obtain 



M s ,s z ,s z = J- Y e m (p)[< S'S' Z \0\1S Z > 5 m0 V2e + T / r 2 dr4 (r)j e (-kr)S(r) 

Jo 2 



L,a 



roc ( / Ax 

y r 2 drD(r)4>f (r)j L ^ — 



or, 



(12) 



+ e i < ^ 2 piia > x: ^-^r 2 '7' o c°; 2 °;°A )L (5')] as) 



with the definitions 



5,(5') = v / 2TTT ^ r 2 dr4 (r)u{-)S(r) 



(14) 
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If we define 



we can write 



Df iL = VWTT r 2 drD(r) ( f ) f{r)j_ 



(15) 
(16) 



M S',S' Z ,S Z = ^ym (p) 



S e 5 mfi < S'S' Z \0\1S Z > +V5J2 C W < S'S' z \0\la > B\ n 



where the explicit dependence on S' of Se and B™ has been suppressed. 
We define the spherical components of the photon polarization vector as 



(17) 



(18) 



The spin structure of the radiative capture operator is assumed to be of the form 

O = <7i • e (19) 

where particle 1 is the proton which is transformed into a A. The expectation values 
of this operator for the baryon spin states are 



Triplet 




< l,0|O|l,l > 


e+ 


< 1,0 \o 1,-1 > 


— e~ 


< 1,0 \o 1,0 > 





< 1,1 \o 1,1 > 


e° 


< 1,-1|C|1,-1 > 


-e° 


< l,l|O|l,0 > 


— e~ 


< l,-l|O|l,0 > 


e+ 


< 1,1|C|1,-1 > 





< 1,-1|C|1,1 > 





Singlet 




< 0,0|O|l,l > 


-e+ 


< 0,0 \o 1,-1 > 


— e~ 


< 0,0 \o 1,0 > 


e° 



The transversality condition gives e° = 0. 

Because of the averaging over the direction of the unobserved momentum p, each 
term in t and m will contribute incoherently. We now consider each of the cases m = 0, 
m = ±1 and m = ±2 separately. In the square-averaging over the photon polarizations 
we omit the uniform factor of 1/2. 
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2.2.1 m=0 

For m = 



oo 

M2' Sz =< S'S' Z \0\1S Z > X>°(P)(^ + V5Cf '° 2 A 1 ^°) (20) 

1=0 

This term will contribute to the transitions (S z = ±1 — > = 0) and (S^ = — ► = 
±1) 

oo 

M^o" 1 = -e- £ *?<P)(Si + -7=^°) = M^ 1 (21) 
1=0 V2 

Mi££ = e + E y Ap)(^ + -^ ) = -<='o (22) 

oo oo 

= -e- £ 3?(p)($ - v^B?); M^'° = e+E l?(p)(S* - v 7 ^ ) (23) 
i=o e=o 

giving a contribution to the total capture rate of 

oo -. 

£(^ + -^°) 2 (24) 
1=0 V2 

from each of the initial magnetic states ±1 for either the singlet or triplet final state 
and 

oo 

2Y / (Se-V2B° e ) 2 (25) 
e=o 

from the initial magnetic state for the triplet final state. 

Since the singlet and triplet states add constructively for S z = 1 and destructively 
for S z = —1, they are incoherent if and only if the populations of the these two initial 
states are equal. 

2.2.2 \m\ = 1 

For \m\ = 1 



l=i 



M^ X " S ' =< S'S' z \0\\S z - 1 > VEcl'-^'f' ]T Yl(p)Bj (26) 
contributes to (S z = — > S' z = 0); 

= -€ _ V?E tf(p)*i = M -=i ( 27 ) 



3 00 
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and {S z = 1 -> S' z = ±1) 

/7y oo rrr oo 

^Y e \ P )Bh M^ = -e+J S -5:Yl( P )B} (28) 
i=\ t=i 

while 

oo 

M^lf' =< S'S' Z \0\1S Z + 1 > V5C£' + £>-l> s * Y f\p) B l (29) 

i=\ 

contributes to (S z = — ► = 0); 

, = e +Vf £ y r i(p)Bi = -M^ (30) 

and (,S Z = -1 5^ = ±1) 

^E^Cp)^ 1 ; ^-^-^^^^(p)^ (3i) 

Z €=1 V Z *=1 

giving a contribution from |m| = 1 to the capture rate of 

=i 



2 

3 H (pi) ' f ° r each ° f Sz = °' ±l ( 32 ) 



Only the S z = state contributes for the singlet final state. Since the relative sign of 
the singlet and triplet is different for m = 1 and m = — 1, the singlet and triplet states 
are always incoherent for \m\ = 1. 

2.2.3 \m\ = 2 

For |m| = 2 

oo 

M^g A =< & S^OllS, - 2 > v^cf ~f 2 f' ]T *7(p)£f (33) 
contributes to (S^ = 1 — > S 1 ^ = 0) 

oo 

M^l = -e-VSY, Y?(p)Bi = M° m % (34) 

1=2 

while 

oo 

M^lf =< S'S' Z \0\1S Z + 2 > y/Ecl' + 2 2 '\ 2 ' Sz Y i~ 2 iP)B 2 e (35) 

1=2 
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contributes to (S z = — 1 — ► S' z = 0). 



£=2 

giving a contribution to the capture rate for either initial spin projection ±1 of 

oo „ 
1=2 



(36) 



(37) 



Since the singlet and triplet states add constructively for S z = 1 and destructively 
for S z = —1, they are incoherent only if the populations of the these two initial states 
are equal. 

2.3 Summary 

Combining the expressions in the previous section we obtain the following results for 
the capture rate for the various initial and final states. The notation for quantities cal- 
culated with the singlet or triplet final state on the A-n system is now shown explicitly. 

S' s 

For the amplitudes in the form Mm' z 

m = 0, e~ 



M, 



o,-i 



-]T yAp) [Si(i) + + s e (0) + -^(0) 

<° = - £ *?(P) (^(1) - V^B?(1)) 



m = 0, 



(38) 
(39) 



Mo^EyAp) 



£=0 



5,(1) + - 5,(0) - -Lfl?(0) 



^o _1 '° = E y Ap) (#(1) - V^B?(1)) 



£=0 



m = 1, e 



M 



0,0 



^/Ie^p)^ 1 ) 



(40) 
(41) 

(42) 
(43) 
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Mr^-vfz^P^Ki) 

m = —1, e~ 



M\- 1 = J-J2Y e - 1 (p)B° i (l) 



m = — 1, e 



M -f = V?E^ _1 i(p) ^(i)-^(o) 



/3 



m = 2, e 



M^ 1 = -V3^y/( P ) [sj(l) + £](0) 



m = —2, e 



+ 



For the capture rates we have 



00 r i i 



1 2 oo 
+ 

£=1 



OO 



+ 3 +*f(0) 



£=2 



S, = 0; 2£[#(l)->/2B?(l)] + 3]T[ J B £ 1 (1)] +3^[^(0) 

£=i i=\ 



s, = -i; E 



S,(1) + -^B?(1) + S,(0) + -Lb?(0) 



2 oo 
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+ 3 -*f(0) 



(52) 



£=2 



Using plane waves for I > 1 



S z = l; 



Sb(l) - Sb(0) + -7= [b °(1) " B§(0)] +3$: + 12 £ [B?] (53) 

V 1 J £=1 £=2 



r n2 r i2 r 

S z = 0; 2 [5 (1) - V2B§(1)] +2^ [#(1) - V2B° e (l)\ + 6j2[ B e 



e=i 



(54) 



£=1 



1 



S z = -l; \S (1) + S (0) + -^[B°{1)+B°{0) \ +4^ 



2 oo 



£=1 



1 



2 oo 



£=1 

(55) 

Adding, and assuming all magnetic projections of the deuteron to be equally populated, 
we find 



s (i) + ^sS(i) 



+ 2 



5 (0) + 



i 2 



+ 2 [S (l) - V2B°(1) 



oo , „n oo „ oo 

l=\ L J £=1 £=2 



Also useful are the expressions for a given photon polarization. 



(56) 



S z = l; 



Eh 1 



£=i 



+ 12£[b? 

£=2 



(57) 



2 °° r n2 °° r 

,S 2 = 0; S (l) - V^B§(1) + £[#(1)-V2B,°(1)] +6E[ 5 1 



e=i 



(58) 



S 2 = -1; 



S (1) + So(0) + ^=[B (1)+B° (0)\ +4^ 



2 oo 



S t + ±B° e 



-,2 o oo 



(59) 
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S z = l; 



S (l) - 5b(0) + ^= [B °(l) " 5 °(0) 



+ 



S 2 = 0; 



2 00 r 

5 (1)-V^B °(1)| +EN 1 )-M°( 1 I 



^ = -1; lE[^ 



(60) 
(61) 
(62) 



2.4 Discussion 

The ratio of the triplet to singlet rates, in the limit of s-wave contributions and an 
s-wave deuteron only, and assuming equal singlet and triplet scattering in the s-wave, 
is two. The introduction of the higher partial waves, and a difference in the scattering 
lengths, modifies this ratio but it is still roughly two. Hence there is a greater sensitivity 
to the triplet scattering length than to the singlet. 

The S and D states of the deuteron interfere coherently, with different coefficients 
depending on whether the initial state is S z = or S z — il. For this reason the shape 
of the spectrum is different for the transition from the S z = sub-state than from the 
S z = ±1 states due to the D-state contribution. Figure [l] shows a typical result for the 
pure I = contribution calculated with an asymptotic A-n wave function. 



3 Solutions with Exponential Potentials 

Rijken et a/.|| have recently fit potential models to the A-nucleon scattering data. We 
use the phase shifts for the s wave determined by this group to define our potentials so 
that the asymptotic form which is singular at the origin can be replaced with a more 
realistic wave function. 

We will represent the effective interaction as a sum of exponential potentials. To 
this end we write the true potential as a Laplace transform 

V(r)= / X(fi)e-^ r dfi. (63) 
Jo 

Consider a calculation of this integral using a Gauss-Laguerre integration scheme with 
a 5 point approximation. The integral has maximum precision for only a single value 
of r, and the value that is chosen as typical for r represents the scale at which the 
integration is made. We choose to work, primarily, at the scale of 1 fm, in which case 
the points in the integration scheme (normally dimensionless) represent inverse ranges 
in units of fm -1 . The smallest of these inverse ranges has a value of 0.26356 fm -1 . 
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Since no exchange with such a small mass is believed to take place in this reaction, 
one expects to find zero for the coefficient of this term. Indeed, fits to the phase shifts 
produce very small numbers. We set these values to zero and consider fits with only 4 
parameters. Thus, the potential is parameterized with the form 

4 

V(r) = V~ wr - (64) 

i=l 

In this case the longest range entering into the problem is that of the second Gauss- 
Laguerre point, which has an inverse range of 1.41340 fm -1 and corresponds very well 
with two pion masses. The other ranges are 3.59642 far 1 , 7.08581 far 1 and 12.64080 
far. We have also considered the scale of one-half fm (which means that the inverse 
ranges are doubled in value) and 2 fm (which means that the inverse ranges have been 
multiplied by 1/2). While we work primarily with the fits on the scale of 1 fm, the 
alternate fits provide estimates of the model dependence. It is the 2 fm range fit which 
is more useful for reasons discussed in section ^. 

3.1 Jost Solutions 

Consider the solution of the Schrodinger equation for a sum of exponential potentials: 

N 

V(r) = X je -^ r . (65) 
i=i 

Jost p3] writes the solution for the s-wave, f(p,r), as 

/(p,r) = e-^£C>)e- m ^ (66) 

7 

where the subscript 7 is a compound index representing a set of N integers. For 
example, for a three-term potential 

7 =[j,M], i,M = 0,l,2.... (67) 
The coefficients C 7 (p) are given by the recursion relation 

, , AiC b -_ liM (p) + \ 2 C^ k - lA {p) + A 3 C b - fc , f -i](p) 
u ' ' 1 m^(m y + zip) 

where 

mj = mij jk fi=jfii + kfi2 + ln3- (69) 
The recursion is started with 
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C[0,0,0] — 1, C[_ ljfcji ] — Cy_ij] — C^jfc^i] — 



(70) 



and is built up by first computing all coefficients with the sum of indices equal to one, 
then two, etc. with no negative index. 

The solution with the proper boundary condition at the origin for an incoming 
spherical wave with unit amplitude at infinity is 



ip(p, r) 



f(p,r) - S(p)f{-p,r) 



where 



S(p) 



2ipr 
/(P,0) 



(71) 



(72) 



These expressions can be used to calculate the values of the S-matrix for any value of 
p. 

In order to calculate the overlap with the deuteron wave function, we require the 
wave function for real (positive) values of p. In this case we can write (for real Xj) 



f(-p,r) = f(p,r) with S(p)=e 2i5 ^. 



We can now write the wave function (Eq. |7l| ) as 

e iS(p) 



tp(p, r) 



2ipr 



-iS(p) —ikr 



e iS(p) e ikrJ2a ( (-p)e- 



(73) 



(74) 



Note that the lowest order term is given by 

pr 



(75) 



which is identical to the asymptotic wave function. We shall use these results to 
represent the A-n wave functions. 



3.2 Fits to the Phase Shifts 

Potentials were fit to the phase shifts(2*|] for each of the six cases of Ref. 0, for both 
the singlet and triplet states and for each of the scales mentioned above. Even though 
the interaction of interest is for A-neutron scattering, the fits were made to the A- 
proton phase shifts since they were more closely related to data. The results are shown 
in Fig. |2] for the triplet case. The potentials themselves are shown in Fig. ||. Table || 
lists the strengths from the fit for the scale of 1 fm. It is this fit which is used in the 
remainder of the paper unless otherwise noted. 



15 



Case 




Ai 




A2 A3 


A 4 


a 


-0. 


.26431 


-45. 


.03399 


304.25626 


2327 


,48169 


b 


-0 


.25513 


-52 


.52803 


352.90689 


3036 


.24731 


c 


-0 


.30427 


-59. 


.00891 


396.75662 


3540 


.90601 


d 


-0 


.26357 


-72 


,61646 


499.84998 


5000 


,94873 


e 


-0. 


.26643 


-79. 


.84088 


561.22314 


5784. 


,05615 


f 


-0 


.27467 


-87 


,65697 


643.43768 


6791 


.61670 










Triplet 






a 


-0 


.43381 


-35 


.94874 


-0.17723 


2502 


,01367 


b 


-0 


.42795 


-37 


.85583 


9.45284 


2803 


.92017 


c 


-0 


.43581 


-39 


.68054 


32.58380 


2944 


,07300 


d 


-0 


.44061 


-41. 


.15301 


52.91377 


3138 


.80127 


e 


-0. 


.40514 


-43. 


.49734 


77.60744 


3344 


.70093 


f 


-0 


.35775 


-46 


.00097 


104.79501 


3541 


,92969 



Table 1: Values of the parameters used in the fit for the scale of 1 fm. 









Singlet 






Case 


a (Rijken et al.) 


a (1 fm) 


ro (1 fm) 


a (2 fm) 


r (2 fm) 


a 


-0.71 


-0.73 


6.71 


-0.63 


4.58 


b 


-0.90 


-0.92 


5.57 


-0.82 


3.99 


c 


-1.20 


-1.23 


4.55 


-1.11 


3.52 


d 


-1.71 


-1.74 


3.71 


-1.59 


2.99 


e 


-2.10 


-2.16 


3.46 


-1.97 


2.57 


f 


-2.51 


-2.59 


3.29 


-2.36 


2.68 








Triplet 






a 


-2.18 


-2.21 


2.99 


-2.07 


2.57 


b 


-2.13 


-2.16 


3.05 


-2.03 


2.71 


c 


-2.08 


-2.10 


3.15 


-1.97 


2.68 


d 


-1.95 


-1.97 


3.33 


-1.84 


2.84 


e 


-1.86 


-1.86 


3.40 


-1.74 


2.94 


f 


-1.75 


-1.74 


3.42 


-1.65 


3.03 



Table 2: Values of the scattering length (a) and effective range (r ) from the fits described. 
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Table || summarizes the scattering lengths and effective ranges corresponding to 
these fits. Also given are the scattering lengths found by Rijken et a/. with the original 
potential. 

In Fig. U are plotted the full wave function and the asymptotic wave function for 
the singlet "a" case of Ref . || . 

4 Results 

We now turn to problems in the analysis and some possible solutions. The first issue 
is the use of the effective range expansion (ERE). In the ir~d — * nwy measurement the 
ERE is adequate to describe the s-wave phase shift over the full range of data being 
analyzed. As is shown in Figure §, this may not be the case for the nA measurement. 
Certainly it is not adequate over the full range, and it may be somewhat questionable 
even if only the upper 10 MeV of the spectrum is used. Of course the ERE is only 
useful if the asymptotic wave functions are used or if the short-range corrections are 
made with a technique similar to that employed in Ref. |§. 

To find an alternative to the ERE we examined the parameters which came from the 
fitting of the exponential potentials to the phase shifts of Rijken et al. It was found that 
an acceptable fit could be obtained with the longest range parameter held fixed. Plots 
of the other three parameters (see Fig. [6| for the triplet case) show a linear behavior 
with case number. We therefore defined continuous variables, which reproduce the 
phase shifts (to a good approximation) at the integers 1 through 6 (corresponding to 
cases a through f), and create an interpolation and extrapolation procedure. 

We define strengths 

Xi = a + xdi (76) 

for i=2, 3, 4. Here x = x s for the singlet case and x = xt for the triplet case. 
The values of Cj, d{ and Ai are given in Table ||. While this method allows a range of 
scattering lengths to be produced with an appropriate short-range wave function, there 
are limitations. Since we would like to be able to study scattering lengths outside the 
range given in Ref. Q, we wish to consider values of x outside the range (1-6). Indeed 
this does give an extension of the range of scattering lengths, but for the 1 fm scale the 
range is limited. Outside this range the scattering lengths return to previous values. 
This constraint dictates the limits of the analysis shown in the figures to follow. It 
was found for the fits on the scale of 1/2 fm that the problem was greater, so that the 
analysis would be restricted to an even smaller area. The 2 fm range fits were used 
for the model dependence estimates because, in this case, a larger range of scattering 
lengths was feasible. However, the scattering lengths do not reproduce those of Ref. 
H as well as the 1 fm scale, as may be seen in Table ^. 

In performing the analysis one may use the asymptotic wave functions or the full 
wave functions. The use of the full wave functions should be preferable, but one can 
obtain an idea of the sensitivity to the details of the full wave function by comparing 
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Al C2 d 2 C3 ^3 C4 ci 4 

Singlet -0.26431 -36.410 -8.624 235.36 68.90 1457.6 869.9 
Triplet -0.43381 -34.576 -1.373 -16.90 16.73 2321.4 180.6 

Table 3: Values of the parameters for the linear representation of the potentials. 



ase 




Ai 




A 2 


A3 


A 4 


1 


-0 


26431 


-45 


03399 


304.25626 


2327.48169 


2 


-0 


26431 


-52 


14749 


349.97018 


3001.88623 


3 





26431 


-61 


00436 


417.78851 


3645.56372 


4 


-0 


26431 


-72 


61623 


499.85406 


5005.63232 


5 





26431 


-79 


96191 


562.18121 


5807.02148 


6 





26431 


-88 


15266 


648.76056 


6839.37939 










Triplet 




1 


-0 


43381 


-35 


94874 


-0.17723 


2502.01367 


2 





43381 


-37 


67183 


9.01196 


2795.74268 


3 


-0 


43381 


-39 


68257 


32.59604 


2939.50488 


4 





43381 


-41 


09288 


52.92378 


3138.51001 


5 


-0 


43381 


-42 


13403 


69.27229 


3286.15967 


6 





43381 


-42 


81454 


83.45535 


3404.96362 



Table 4: Values of the parameters obtained for the 1 fm fit with constant \\. 



results with those for the asymptotic wave functions. Figure [7| shows the sum of the 
S z = and S z = ±1 s-wave contributions to the spectra as shown in Fig. [I] for the 
asymptotic wave functions (solid lines) and the full wave functions (dashed lines) for 
the parameters of the singlet "a" and "f" cases. 

Perhaps the largest problem in the measurement is the separation of the singlet and 
triplet scattering lengths in the final state. If no spin degrees of freedom are measured 
then the singlet and triplet states contribute incoherently to the rate, and it might 
seem impossible to separate them. However, because the interference between the S 
and D states of the deuteron differs, the shapes of the singlet and triplet spectra are 
different, as illustrated in Fig. ||. Because the triplet gives a larger contribution than 
the singlet, a greater sensitivity to the triplet scattering is seen in the fits of the data 
and pseudo-data in the figures which follow. 

Figure |9] shows an analysis of the data of Gall et al. [16]. Since this experiment 
was only a feasibility study, one can not expect to obtain much information about the 
scattering lengths. It is interesting to see, however, that the values are in the range 
expected. The open circle in the upper left hand corner shows the point of minimum 



X 2 



Figure |l^ shows a similar analysis of the Gall et al. data with the 2 fm scale fits. In 
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this case a larger range of scattering lengths can be studied so that an error estimate 
can be given. One-standard-deviations limits can be read from the graph (the inner 
contour) with the triplet length lying between -1.3 and -2.6 fm and the singlet length 
between -0.2 and -6.3 fm. 

Figure 11 provides a comparison of the experimental and theoretical spectra for 
the scattering lengths corresponding to the minimum x 2 (=37.8 for 37 data points) at 
a s = —2.96 and at = —1.72 for the 1 fm scale. 

Figure 12 presents a similar analysis of a pseudo data spectrum. The pseudo data 
were generated from a selected spectrum by including errors chosen from a Gaussian 
distribution such that they have a value of 3% of the rate at the maximum and are 
proportional to the square root of the rate at other points, as would be the case with 
errors dominated by counting statistics. The analyses shown were made over the full 
range of the Gall et al. data from 255 to 293 MeV unless otherwise specified. The case 
chosen for presentation in Fig. [l2| is a favorable one in the sense that the scattering 
lengths used for generation of the pseudo-data (solid circle) fall near the center of the 
inner ellipse. As is necessary statistically, most of the results do fall nearer the edge of 
this ellipse with about 1/3 falling outside. One can observe from this figure that the 
experimental uncertainty in such an experiment would be ~ ±0.3 fm for the triplet 
scattering length and ~ ±0.8 fm for the singlet. 

It is possible to use such an analysis of pseudo data to estimate errors due to model 
assumptions. For the following we use the same pseudo data that were described 
above. Figure 13a shows the results of an analysis using asymptotic wave functions 
with phase shifts determined from the Jost solutions. Figure 13 b shows an analysis 
with the additional approximation that the phase shifts are taken from the ERE. As 
can be seen a significant error results. 

Use of the asymptotic wave functions is an extreme approximation. Although we 
do not know exactly which more realistic wave function should be used, we can choose 
one among the Jost models we have been considering. To this end we analyze the 
spectrum created with the 1 fm scale using the 2 fm scale spectra. The result can 



be seen in Fig. 14. The positions of the singlet and triplet scattering lengths at the 



minimum are -1.613 and -1.884 fm to be compared with the values found with the 
(1 fm scale) analysis above of -1.896 and -2.031 fm (the input values are -1.968 and 
-1.827 fm). Thus we find model errors of ±0.28 fm for the singlet and ±0.15 fm for 
the triplet case. We will take these as estimates of the model (theoretical) errors. 

Figure 15 shows the same type of analysis but including only the upper 10 MeV of 
the spectrum in the comparison. It is seen that the results are very similar, due perhaps 
to the fact that the main difference in shapes between singlet and triplet occurs in this 
region (see Fig. ||). The analysis over this range may offer experimental advantages in 
addition, since the background due to neutral pion decay is less severe. 

We have seen that it is difficult to determine the scattering lengths individually. It 
is possible to separate the singlet and triplet states with the use of spin information 
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from the deuteron initial state or the photon final state. Figure ^ shows the spectra for 
the six possible combinations of these spin projections for four selected parameter sets. 
While the overall normalization of these curves is arbitrary, the relative normalization 
among them is correct. The upper left hand curve (for S z = +1 and right circularly 
polarized photon) shows the maximum sensitivity since it corresponds to the case 
where the singlet and triplet final states are coherent in a destructive manner. The 
disadvantage is that this same cancellation leads to a small rate if the singlet and 
triplet scatterings are not too dissimilar, as can be expected from SU(6) symmetry. 
However, this small quantity is a direct measure of the difference between the singlet 
and triplet lengths. Note that the spectra from an initial state S z = are independent 
of the singlet parameters. 

We can express the results in terms of more conventional variables. Figure |l~7| 
illustrates the spectrum expected from a pure vector polarized deuteron target with 
the asymmetry defined as 



T(S Z = +1)-T(S Z = -1) 
T(S z = +l) + T(S z = -l) 



We see that a measurement with precision of the order of 2% is needed near the peak 
of the spectrum or 10% at the upper part of the spectrum. What is being searched 
for is a deviation of the asymmetry from -1 or, in other words, a non-null value of the 
rate from the initial state S z = +1. Such an observation would prove that the two 
scattering lengths are unequal. The dashed curve corresponds to a difference of 0.1 fm, 
the dash-dot curve to 0.66 fm, the solid curve to 0.91 fm and the dotted curve to 1.46 
fm. This result is almost completely free from the model dependent effects discussed 
previously. The observation of a difference from -1 in the asymmetry would be a direct 
indication of a difference in the scattering lengths. 



Figure 15 shows the spectra expected from measurements in which the circular 



polarization of the final photon is measured. The asymmetry is defined as 

_ r( £ = +i)-r( e = -i) 
^ - r(e = +i) + r( e = -i)" (78) 

Here again the same order of accuracy is needed to separate the dotted and dashed 
curves. 



5 Conclusions 

A formalism has been presented with attention to the spin degrees of freedom. Due to 
the difference in shape of the spectrum for the reaction proceeding from the zero mag- 
netic quantum number projection of the deuteron from that with magnetic quantum 
numbers of ±1, the singlet-triplet separation can be made, although the difference is 
not large and the triplet state tends to dominate. 
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An analysis of the present data which exist for this reaction [16] has been made. 
Even though the uncertainties in the data from this feasibility experiment are large, 
some information can be obtained; in particular, separate numbers can be extracted 
for the singlet and triplet scattering lengths (with large overlapping errors). We find 
for the singlet scattering length in the range —0.15 to —5.0 and the triplet value from 
— 1.3 to —2.65 (from the 2 fm range fit). 

While this precision may seem rather modest, it is useful to compare with current 
values in the literature. Tan[24| obtained a value (believed to be mostly the triplet 
scattering length) of —2.0 ± 0.5 fm from an analysis of data on the reaction K~d — > 
■K~pA. Recent data|l7] on the reaction pp — > pK + A, in conjunction with previous 
elastic scattering data[|], lead to a value for the spin-averaged scattering length of 
—2.0 ± 0.2 fm[17|. The importance of the influence of the third strongly interacting 
particle in the final state is difficult to estimate and the errors are experimental only. 

In a three-body calculation for the hyper-triton bound state Miyagawa concludes |25[| 
that the best values of the singlet scattering length lie between —2.7 fm and —2.4 fm 
and the triplet between —1.6 fm and —1.3 fm, conforming to the "f" solution of Rijken 
et aZ.gj. 

An analysis of pseudo data for a K~d — > Awy experiment with a cross-comparison 
of models leads to an estimate of the model dependence of the order of ±0.2 fm. 
The same type of study gives an estimate of the error in such an experiment of ±0.3 
fm for the triplet scattering length and ±0.8 fm for the singlet. These estimates are 
for measurements without spin information and assume an uncertainty of 3% for the 
maximum rate in the spectrum. 

It was shown that spin information in either the initial or final state would be 
valuable in separating the scattering lengths. If it were possible to perform the capture 
experiment on a deuteron target of purely magnetic quantum number zero the triplet 
scattering length alone would be measured. 

The use of a polarized deuteron target would also allow the separation since the sin- 
glet and triplet states interfere destructively or constructively according to the relative 
alignment of the deuteron spin along or against the direction of the photon. 

Since the previous feasibility study showed that the measurement was possible 
but very difficult, one might ask if the increased difficulty due to the requirement 
of polarizing the deuteron might not make it impossible. Certainly it adds another 
constraint, especially if it is necessary to add high-Z material which would preferentially 
capture the Kaons. In this regard we mention a possible polarized deuteron target 



which is made from a hydrogen-deuteron molecular system [26, K7j. It has a high 
density and contains no heavy materials. 

For a measurement of the circular polarization of the photon, the spectra show a 
significant sensitivity of a similar type. 

This work was supported by the U. S. Department of Energy and the National 
Science Foundation. HKH thanks the Korean Scientific and Engineering Foundation 
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- a =— 2.2 fm, r =3.2 fm - 
1.2 r ° - 

1.0 1 / \ - 




255 265 275 285 295 

k (MeV/c) 

Figure 1: Comparison of the spectrum from the S z — state with that from the non-zero 
initial projections. The dashed curve is the spectrum from the S z = magnetic sub-state of 
the deuteron. 
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Figure 2: Comparison of the fits for the s-wave phase shifts with the results of Ref. || for 
the triplet states. The results for the scale of 1 fm is represented by the solid line and the 
scale of 2 fm by the dashed line. 
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Figure 3: Potentials obtained from the fits for the triplet case. The meaning of the solid and 
dashed lines is the same as in Fig. |2|. 
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Figure 4: Wave function for the asymptotic wave function (dashed curve) compared with 
the Jost solution (solid curve). 
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Figure 5: Phase shifts from the fitted potentials (solid) and from the effective range expansion 
(dashed) for the case c of Rijken et al. for the singlet and triplet cases. 
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Figure 6: Parameters from the fit to the phases shifts of Rijken et al. showing the linear 
representation for the triplet case. 
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Figure 7: Comparison of the spectra expected (£ — contribution only) with an asymptotic 
wave function (solid curves) with that using a solution to the Schrodinger equation (dotted 
curves) for the Jost parameters corresponding to cases "a" and "f" for the singlet final state. 
The lower curves are those resulting for the singlet spectrum and the upper ones are for the 
triplet spectrum. 
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Figure 8: Comparison of singlet (dashed) and triplet (solid) spectra for a full calculation 
with realistic wave functions. The dash-dot curve shows the singlet spectrum renormalized 
to the triplet at the lowest energy. 
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Figure 9: Analysis of the data of Gall et al. using the 1 fm scale fits. The contours are lines 
of constant y 2 . Points lying inside the curve labeled +1 are within one standard deviation 
of the minimum. The solid squares are the values taken for the 6 cases of Ref. || and the 
open circle is the value at the minimum \ 2 of the fit. 
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Figure 10: Analysis of the data of Gall et al. with the 2 fm range fits. The innermost curve 
corresponds to one standard deviation, the next to two, etc. 
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Figure 11: Comparison with the data of Gall et al. with the best fit using the 1 fm scale. 
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Figure 12: Analysis with a full Jost model of a pseudo data set generated with 3% errors. 
The 1 fm scale fit was used to generate the data and make the analysis. The solid circle is 
the input value chosen for the test and the open circle is the central value corresponding to 
the minimum \ 2 ■ 
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Figure 13: Analysis of the same pseudo data set as in Figure [12] with asymptotic wave 
functions. In the upper part the correct phase shifts are used while in the lower portion the 
effective range expansion is used. 
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Figure 16: Spectra for the various spin combinations of the initial deuteron and the final 
photon. The three figures to the left are for right circular polarization and those to the 
left for left circular polarization of the final photon. The three rows are for S z = +1, 
and —1 from the top down. The solid curve has a s = —1.15, a t = —2.06 the dashed curve 
a s = —1.96, a t = —2.06, the dash-dot curve a s = —1.15, a t = —1.81 and the dotted curve 
a s = -0.60, at = -2.06 fm. 
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Figure 17: Vector polarization asymmetry for the deuteron initial spins. 
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Figure 18: Circular polarization of the photon under the same conditions as the previous 
figures. 
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